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Abstract The open-flavor strong decays are studied based on the interaction of potential quark model. 
The decay process is related to the s-channel contribution of the same scalar confinement and one-gluon- 
exchange(OGE) interaction in the quark model. After we adopt the prescription of massive gluons in time-like 
region from the lattice calculation, the approximation of four-fermion interaction is applied. The numerical 
calculation is performed to the meson decays in u, d, s light flavor sector. The analysis of the D/ S ratios of 
bi — )■ LOTi and ai — >■ pn shows that the scalar interaction should be dominant in the open-flavor decays. 
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1 Introduction 



Although QCD is considered a correct theory 
for strong interactions, knowledge about the hadron 
structure in low energy region is restricted due to 
the color confinement. The potential quark model is 
widely used to identify the conventional hadron states 
in hadron physics. In contrary to its impressive suc- 
cess in hadron spectra, especially when heavy quarks 
are involved, its interpretation to the hadron strong 
decays is unsatisfied. Although there exist some phe- 
nomenological models for the strong decays, the re- 
lationship between these models and the potential 
quark model is somewhat obscure. 

One of the most popular models for open-flavor 
strong decays is ^Pq model developed in the 1970s 
[l|, 0] ■ This model has successfully demonstrated its 
universal practical utility when applied to a great 
number of particular decay channels Later 
the flux-tube-breaking model was proposed and ^Pq 
model could be regarded as a limiting case of this im- 
Toved model As early as 1978, Eichten et al. 
developed the Cornell model by incorporating 
the possibility of creation of a light-quark pair into 
the quark model Hamiltonian. However, in their 
model they considered the quark interaction the time- 
component part of the vector interaction and assumed 
that the interaction of the quark pair creation was 
the same as the instantaneous interaction between 



two constituent quarks. In recent years an extended 
model including the scalar confining and vector OGE 
interactions was studied by E.S.Ackleh et al. 10|. A 



similiar scalar color-singlet confining interaction was 
dervied from a relativistic string breaking kernel ll| . 

The instantaneous interaction in the above mod- 
els always assumes Breit approximation when deal- 
ing with gluon's momentum. For the potentials in 
quark model, the energy of the exchanged gluon is 
negligible as compared with the the masses of con- 
stituent quarks. Therefore the transferred gluon mo- 
mentum is space-like. Nevertheless, it is in all the 
probability time-like if considering the creation of 



Besides, based 



quark-antiquark pair by the gluons. 
on the recent study in lattice field theory 12|,ll3|, glu- 
ons are supposed to act as massive vector bosons in 
non-perturbative region with masses evaluated about 
600^1000MeV. A non-vanishing gluon mass is also 
needed in the phenomenological calculation of the 
diffractive scattering IJ] and radiative decays of the 
J/V' and T [lil. 

In this paper, an alternative study of the open- 
flavor strong decays is taken and inspected by the 
experimental decay widths. Following Ref. [10|, the 
quark pair-creation interaction consists of a scalar 
confining interaction and an OGE part. We will dis- 
tinguish the Breit approximation of the gluon's prop- 
agator in the time-like region from that in the space- 
like region. In the time-like non-perturbative region, 
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the massive gluon prescription is adopted according 
to Refs. 12, Il3|. In this way, the decay interaction 
will be further simplified to the form of four-fermion 
interaction. 

2 The Decay Model 

To describe the creation of a light-quark pair in 
the quark model, a plausible approach is to consider 
the field quantization of the quark potential. In the 
Cornell model, the quark potential is replaced by an 
instantaneous interaction 1 91 



where 



(i'x(fy:pa[x)-V{x-y)p,{y):. (1) 



(2) 



is the quark color-charge-density operator. Here ip{x) 
denotes the quark field with fiavor and color indices 
suppressed, and stands for the Gell-Mann ma- 
trices for SU{3) generators. Since the confinement 
should be the Lorentz scalar, in Ref. [l^ the instan- 
taneous interaction is replaced by the combination 
of the scalar confinement interaction and the vector 
OGE interaction. 

We will start from the covariant nonlocal current- 
current action of the quark interaction [l6j |: 



d'^xd'^y%l){x)^^T^il){x)G{x-y)'tj){y)^^Tai){y)^ 



d'^xd'^y'4){x)Taip{x)S{x-y)i){y)Tai){y). (3) 



The vector kernel G{x — y) corresponds to the 
gluon propagator in coordinate space which gener- 
ates the OGE Coulomb potential — in the quark 
model. In the momentum space 

Anas 



(4) 



On the other hand, the scalar kernel S{x — y) should 
generate the linear confining potential |6r. Thus in 
the momentum space 

Si,^) = -'-^- (5) 

The relevant coupling constants and b are the po- 
tential parameters in the potential quark model. 

The lattice calculation shows that the behavior 
of the gluon propagator is quite_different in the non- 
perturbative region. In Refs. 
propagator is assumed to be: 



12 



13|, the transverse 



Z{q') 



(6) 



q-~AP{q^)' 

where M{(f) is the running gluon mass. Then the 
kernels G and S are modified to 



G{q') 



q''-A'P{q^)' 
67r5 



(7) 
(8) 



The lattice simulations suggest M(0) = 600 1000 
MeV which means that gluon gets a non-vanishing 
mass Mg in the non-perturbative region q <C Aqcd ■ If 
the q^ term in the gluon's propagator is neglected in 
the quark-antiquark pair-creation process, then the 
decay interaction is simplified to the four-fermion in- 
teraction. The interaction Hamiltonian density for 
na.ir-c.reatinn turns to be: 



where Tis (x) and H„ {x) represent the scalar and vec- 
tor interaction respectively: 



Hsix) ^^J{x)T,i:{x)i,{x)T'^i^{x) 



(10) 



3 



4>ixh^T^ibix)^{x)j^Tyj{x). (11) 



From the interaction in Eqs. ([T0|) and (fTTj) . we 
can derive the formulae for decay rates within non- 
relativistic limit. As meson states are normalized to 
2E in our work, 



{p\p')^2ES'{p-p'). 



(12) 



the differential decay width in two-body decay pro- 
cess A ^ B + C is expressed in terms of transition 
amplitude as: 



S\M\' 



i2nyd\PA-PB-Pc 



d^PBd'Pc 



where S is the symmetric factor 



l + S{B,C) 



(13) 



(14) 



The amplitude A4 is related to the decay interaction 
through: 



M^<BC\Hi{0)\A>^M,. + Ms 



(15) 



where and A4s are the amplitudes from vector 
and scalar interaction respectively. 

For each interaction, the transition amplitude 
comes from four diagrams. For the vector interac- 
tion, its Feynman diagrams are shown in Figure [TJ 
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Fig. 1. Contributions to vector interaction. 




(a) 




(b) 




(c) 




(d) 



In Diagram (a) the qq pair is created from gluons 
emitted by the initial quark while in Diagram (b) the 
qq pair is created from gluons emitted by the initial 
anti-quark. Diagrams (c) and (d) come from the in- 
terchange of final particles B and C in Diagrams (a) 
and (b) respectively. 

The total decay width is expressed as 



Ml 



LS 



(16) 



where Pf — Pb = —Pc in the rest frame of initial 
particle A and TW^'^ are the partial wave amplitudes. 

For further simplification, the space wave func- 
tions of all meson states are taken to be the sim- 
ple harmonic oscillator (SHO) wave functions with 
a common oscillator parameter /3. The partial wave 
amplitudes M.^^ are presented in the appendix. 

3 Results and Analysis 



In the numerical calculation, the common oscilla- 
tor parameter (3 = 400MeV is adopted from Ref. 16 1 
All related masses of mesons are taken from Ref. 17 1 



Other parameters, like constituent quark masses, cou- 
pling constants and h are also taken from Ref. [l6| . 
They are m„ = = 220MeV, m, = 419MeV, 



0^ = 0.60, and fe = 0.18GeV^ respectively. 

The only one parameter which cannot be deter- 
mined from the quark model is the effective gluon 
mass Mg. In this work, this parameter is fixed in 
a least square fit to the experimental decay widths. 
We find that Mg = 668MeV which falls within the 
range 600 ^ lOOOMeV estimated in the lattice calcu- 
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lation 

The results of the decay widths are tabulated in 
Table [T] together with the decay widths of ^Pq model 
and the experimental data. In the table, Fi indicates 
the decay rates of our calculation. As can been seen, 
the widths of the decay processes characterized by the 
creation of ss pair are rather small compared with the 
experimental data. The reason is that the creation of 
ss pair is suppressed in the four-fermion interaction 
due to the heavier mass of s quark. Note that in the 
^Pq model, the transition operator is independent of 
the flavor mass in qq pair creation. This shows that 
the qq pair creation has the flavor 5't/(3) symmetry. 
In the third column F2 in Table [1] the decay widths 
related to ss pair creation are recalculated with the 
flavor symmetry restored with vris = — 220MeV 
in the decay Hamiltonian. The corresponding decay 
widths are enhanced which improves the fit to the 
experimental data. 



Table 1. The decay widths. The decay widths of ^Po model are taken from Ref. [^. The experimental data 
are taken from Ref. Unit: MeV 
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Channel 


Ti 


r2 


Ta 




Exp 


p — V niT 


109.7 




138 


96 


149 


bl — ^ LUTT 


57.7 




160 


176 


142 


a2 — >■ pTT 


51.5 




49.7 


65 


75.4 


a2 — >■ »77r 


15.4 




19.5 




15.5 




1.64 


9.03 


2.26 


11 


5.24 


0.2 


1.20 




1.52 




0.567 


7T2 — ^ /2 7r 


58.6 




77.1 


147 


146 


7V2 — ^ pTT 


58.3 




128 


232 


80.3 


TT2 — > K* K + C.C. 


0.23 


10.8 


4.66 


38 


10.9 


1X2 — >■ pu) 


9.23 




25.1 




7.00 


p3 — ^ TTTT 


47.3 




59.6 


116 


38.0 


P3 — > WTT 


17.8 




17.2 


36 


25.8 




0.68 


7.06 


0.94 


9.2 


2.54 




136 




172 


203 


157 




1.07 


5.98 


1.48 


7.2 


8.51 


/4 — ^ ojo; 


23.9 




14.9 


53 


54 


/4 — TTTT 


27.6 




34.8 


123 


40.3 




0.18 


3.45 


0.25 


5.4 


1.61 


/o(1500)-i>7r7r 


108 




34.7 




38.0 


/o{1500) ^ii'A' 


4.95 


7.99 


0.49 




9.38 




1.96 


2.63 


2.18 


2.37 


2.10 




84.3 


78.9 


21.0 


117 


64.8 


K* -^Ktt 


41.7 


45.9 


46.4 


36 


50.8 


K* (1410)^ Kit 


21.7 


32.4 


1.15 




15.3 




348 


1062 


194 


163 


251 




81.7 


73 


90.2 


108 


49.2 




23.0 


20.3 


20.6 


27 


24.3 


Ki-^Kp 


7.22 


6.20 


6.42 


9.3 


8.57 




2.12 


1.82 


1.88 


2.6 


2.86 


K*^Kp 


14.4 


9.54 


12.5 


24 


49.3 




19.5 


13.8 


17.7 


33 


31.8 


K'^Kn 


45.1 


32 


49.4 


87 


30.0 




18.4 


10.1 


20 


55 


19.6 




0.37 


2.44 


0.23 


3.2 


2.8 



Individual decay amplitudes from scalar and vec- 
tor interactions are listed in Table [51 The scalar in- 
teraction is dominant in most of the decay channels. 
However in the channels ID^IP+IS, 2S^1S+1S 
and 2P — > 15 + 15 the contribution from the vec- 
tor interaction is important, while in the channel 



^Pq — > ^5o + ^5o the vector interaction becomes dom- 
inant, as can be seen in the process /o(1370) — > tttt 
whose decay width may amount to 1000 MeV which 
is too large compared with the experimental result of 
200-500MeV. 



Table 2. The individual amplitudes from scalar and vector interactions. Unit: MeV 
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channel 


L 


S 


Ms 




M 


p — >■ TTTT 


P 





-1.39 


-0.55 


-1.94 


bl — > tJTT 


S 


1 


2.34 


-0.35 


2.00 


61 — > U)TT 


D 


1 


0.62 


0.49 


1.11 


0.2 pi" 


D 


1 


-1.32 


-0.79 


-2.11 


a2 — rjiT 


D 





-0.73 


-0.29 


-1.02 


a2^KK 


D 





-0.28 


-0.09 


-0.36 


0.2 — >■ '?'•"" 


D 





-0.28 


-0.11 


-0.39 


7r2 -> /2 TT 


S 


2 


2.32 


-5.52 


-3.20 


7r2 — 5- /2 TT 


D 


2 


0.383 


-0.003 


0.380 


7r2 /2 TT 


G 


2 


0.01 


0.004 


0.014 


7r2 — ^ pTT 


P 


1 


1.85 


-0.72 


1.14 


772 ~^ P^T 


F 


1 


1.10 


0.88 


1.98 


■7T2 — > K'K + C.C. 


P 


1 


0.49 


-0.56 


-0.07 


■K2 — > K* K + C.C. 


F 


1 


0.09 


0.07 


0.16 


1X2 — ptt) 


P 


1 


1.39 


-0.09 


1.31 


7r2 — >■ ptD 


F 


1 


0.13 


0.10 


0.23 


P3 — ^ TTTT 


F 





-1.31 


-0.52 


-1.83 


P3 — CJTT 


F 


1 


-0.80 


-0.47 


-1.27 


ps-^-KA" 


F 





-0.18 


-0.06 


-0.24 


/2 — >■ TTTT 


D 





-1.95 


-0.77 


-2.72 


f2^Kk 


D 





-0.23 


-0.08 


-0.31 


/4 — ^ CdCJ 


G 





-0.08 


-0.03 


-0.11 


/4 — ^ CdCJ 


D 


2 


0.88 


0.87 


1.75 


/4 — > wo; 


G 


2 


0.15 


0.15 


1.30 


/4 — TTTT 


G 





-1.09 


-0.43 


-1.52 


/4^i<'i? 


G 





-0.10 


-0.03 


-0.13 


/o(1500)->7r7r 


S 





0.95 


-3.56 


-2.61 




S 





0.13 


-0.77 


-0.64 


4>^K+K- 


P 





-0.40 


-0.20 


-0.60 




D 





-0.84 


-0.28 


-1.12 


K* ^Ktt 


P 





-1.04 


-0.50 


-1.54 


K*{IA10)^K-K 


P 





-0.18 


1.39 


1.21 


K* K-K 


S 





-2.32 


7.19 


4.87 


2 


D 





-1.58 


-0.78 


-2.36 


K*^K*-K 

I 


D 


1 


0.92 


0.61 


1.52 




D 


1 


-0.59 


-0.39 


-0.98 




D 


1 


-0.32 


-0.21 


-0.53 




F 


1 


-0.74 


-0.50 


-1.24 




F 


1 


0.74 


0.44 


1.18 


KI^Ktt 


F 





-1.27 


-0.63 


-1.90 


Kl^K-K 


G 





-0.86 


-0.43 


-1.29 




G 





-0.003 


-0.002 


-0.005 




D 


2 


0.15 


0.15 


0.30 


Kl^K*(t) 


G 


2 


0.007 


0.007 


0.014 


One of the important criteria for the strong de- 1 
cay models is the D/ S amphtude ratios in the decays 
61 — > WK and fli — >■ pTT. Experimentahy, these ratios 


D2 - 


^/2p){bh + 2aml) 

2h{hp) ~ 24/32) _^ 4q,j„2 (p2 _^ 





are = 0.277±0.027 and D2 = -0.062±0.02 respec- According to the preceding values of the model 

tively [l7|. In the current model, analytic expressions parameters, the ratios' numerical resuhs are: = 
for these ratios are: 0.566 and I?2 = 0.731. It is apparent that is about 



two times larger than the experimental value while 

•\/2j>f (156 + 8m^a) the calculated value of D2 has a wrong sign. Since 
(17) 

Ihbp^j _(_ 8m^p'ja — 726/3^ ' ^ ' the /3 value is dependent on the wave function of me- 
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son, in Figure [2] we show the dependence of the ratios 
on (3. 





Fig. 2. DjS ratios for Di and Da- 

With respect to ratio, it decreases with an in- 
creasing /3. When /3 rises up to 524MeV, the ratio re- 
generates the experimental value 0.277. Nonetheless, 
as to D2 ratio, the numerical value keeps its opposite 
sign since this ratio changes rather slowly with /3. 

Based on the fact of the dominance of the scalar 
interaction, a scalar-kernel-scalar i^sKs) decay model 
was proposed As a result we consider only the 
contribution from scalar interaction while leaving the 
vector interaction aside. The best-fitted value for Mg 
now becomes 597MeV and the fitted decay widths are 
listed in the Fa column in Table [TJ 

One of the advantages of considering scalar inter- 
action alone is that the DjS ratios are greatly im- 
proved. Now these ratios turn to be: 

5^/2p2 



y2(5j9^- 24/32)' 



(20) 



We obtain Di = 0.264 and = -0.140 when 
/3 — 400MeV, well fit to the experimental results. 
Another improvement lies in the specific channel, 
^S'o + ^S'o. As in the process /o(1370) tttt, 
with the negligence of vector interaction, the decay 
width becomes a reasonable value 318.5MeV. 

4 Summary 

To summarize, we have studied a decay model 
based on the potential quark model. The model in- 
corporates the decay interactions of scalar and vector 
quark currents which are in accordance with the con- 
fining and OGE potentials in the quark model. In 
the non-relativistic limit, the massive gluon propaga- 
tor is assumed and the decay interactions are reduced 
to four fermion interactions. In this framework, we 
have calculated 34 decay channels. The results fit the 
experimental data comparable to the popular ^Pq de- 
cay model if the S'C/(3) flavor symmetry is assumed 
in the decay processes. Meanwhile the results also 
show the dominance of the scalar interaction in most 
of the decay channels. Besides, the scalar interaction 
is also preferred by the DjS ratios of 61 — >■ w -I- tt and 
ai — > /9 + TT. Thus we have calculated decay widths 
with only scalar interaction (the sKs model which is 
quite similar to the ^Pq model since the qq scalar cur- 
rent produces a ^Pq quark pair). It seems that scalar 
interaction alone is able to give a crude estimation of 
most decay widths. 

We would like to thank Prof. Shi-Lin Zhu for use- 
ful discussions. 



'24/32 



(19) 



Appendix 



where C/ is the flavor factor: 



A Amplitudes and Overlapping Inte- 
grals for Some Channels 

The partial wave decay amplitude A-^ B-\-C can be 
expressed as 



H.LS V8EaEbEc „ ,,LS 



(Al) 



C}^{2Tb + 1){2Tc + 1) 



Ta 



Tb 

t2 



Tc 
ti 



(A2) 



where Ta, Tb, Tc are the iso-spins of mesons A, B, 
C, respectively, ti, t2, t are the iso-spins of quarks la- 
beled as 1, 2, 3 in Figure[T] respectively. Similarly, in the 
following, the masses of quarks 1, 2, 3 will be denoted by 
mi, 7712, n^, respectively. 
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The decay amplitudes can be split into two parts 
which include scalar M^^ and vector M^^ . Furthermore, 
each MI"^ would part into four components: 

M/'^ = M/-^ (a) + Aff (b) + Mi'^{c) + Mi^{d) . (A3) 

according to Figure [T] In the following, we will only 
present the formulae for Mt^{a). Mt^ {h) is related to 
Mi^{a) by a charge conjugate: 

{b' B ^C) — (^^ly^B+Jc-s+SA+Ss+Sc+i 
xMi^{a-A^C + B). (A4) 



where A, B, C are the charge conjugates of A, B, C, re- 
spectively. The Mt^{c), Mt^{d) are related to Mt^{a), 
Mi^ih) by the exchange of final particles B and C. 

Since the decay interaction is four-fermion interaction, 
the spatial overlap integrals involve: 



PAc{mA,mc) = 
VAc{rnA,mc,m) = 



vb (mB,m) 

where 



(fc) 
(fe) 



,(k)km. (A5) 



and 



km = {k^ m = (A6) 

^+^{kx — iky) m = — 1 

All the spatial wave functions tpnim are taken to be the 
simple harmonic oscillator (SHO) wave functions with the 
oscillator parameters /3a, /3b, /3c for mesons A, B, C re- 
spectively. We have 



PBims) ■ 



_ 5{Ib,0) 



(A7) 



Vb {mB,m) = 



S{lB,l)S{mB,m) 3 



ubI 



Let 



(A8) 



(A9) 



Below we list the non- vanishing integrals pac and vac 
relevant to our work. 



PAC (0,0) = r) 
WAc(0,0,0) = - 



-Pf 



• 25^15' 

PAC(0,0) 

«Ac(0,0,0) 

• IP -^15' 

PAc(0,0) 

?;ac(0,0,0) 

UAC(1,0,-1) 

• ID ^15' 

PAC (0,0) 

7;ac (0,0,0) 



(3p\-3l3*c-'2p}Pie)v 
-4/3i/3S + 2p^C'/3i) 



V2riPAPl{P\ + Pl~p}e) 

~V2-nPAPl 



i;ac(-1,0,1) 



2p'//3ig^ 



2pffilPUv{W\ + Wl-p)e 



?;ac(1,0,-1) = «ac(-1,0,1) = 



-2pfPlplr,i 



• ID^IP 



PAC(0,0) 

Pac(1,1) 
«Ac(0,0,0) 
^'Ac(2,l,-1) 

^^AC(1,1,0) 



PAci-l,-l) = 
2V2l3ll3cV 



[2/3S - ^pjPU' +ppU' + f3\m 



«ac(-1,0,1) = 
z;ac(-1,-1,0): 



2./2p\Plr,{P\+Pl-ep)) 

-VAC{IA-1) 
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Vac{0,1,1) = Vac{0,-1,-1) = 

Pac(0,0) 
Wao (0,0,0) 



VAC{-1AI) = 



We further introduce some useful combinations: 



AL=PfPB{0) ^ pAc{mA,mc){lcmcLO\lAmA) 

mAmc 

Blj=^/21a + Ipb{Q) ^ VAc{rnA,rnc,m){lcrncLO\Jmc) {lAmAlrnlJmc) 



The relevant partial wave amplitudes are given in subsections. 
A.l S^S + S 



Mj°(a) 



'V2M^ 
\/27Ta, 



3M^ 



m + TTii m + m,2 m 



mi(mi+m2) -m(mi -Sma) ^^ _^,_3 1_ 

mi(m-|-mi)(m-|-m2) mi m 



Ml\a) = V2Ml°{a;'S^^'So + 'So) 

HI,,,, 2V^a 



/ mi+ma 1 1 

( — ;^ — ; + ( j^n 

mi[m + m2) mi m 



M,"(a) 



Ml\a) = -V2Ml\a-^Si^^So + ^So) 
Ml\a) = -s/2Ml\a-;'S,^^S, + 'S^) 



\/3Mf (a;35i^i5o + i5o 

27ra, 



M„"(a) 



m(3mi-m2) + mi(mi-l-m2) ,1 1 , 

7 — ; TT — ; ^ ^0 - ( 1 j-t>ii 

mi(m-|-mi)(m-|-m2j m mi 



-V3M^»(a;^5i^^5'o + '5o) 



(AlO) 
(All) 
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A.2 P^S + S 



M°°{a) 



V2M^ 

2'Kas 
3Mf 



: + — : ^1 + — Soo 

m + rrii m + m2 m 



mi (mi — m) + m2 (mi + 3m) ^ ,3 ''^ 
mi(m-|-mi)(m-|-m2) ^ mi m °° 



3^/5M2 



m + mi 171 + 1712 ^™ 



2mi(mi-m)-|-2m2(mi-|-3m) /- 3 1 

; r- ; r Vli-|-Vi( )±«22 

mi(m-|-mi)(m-|-m2) mi m 



Mf{a) = ^M'f{a-'P2^'So + 'So) 



Ml\a) = - 



30M2 



2(mi+m2) - , A^. 1 1 ^ o 

; 7-Ai-|-V2( )S22 

mi(m-|-m2) mi m 



• ^1^0 + 3^1 



Mf{a) = -^M^°{a;'P2^'So + 'So) 



Ml\a) 



; ; Ai-|-V2( j±;22 

mi(m-|-mi)(m-|-m2) mi m 



A.3 D^S + S 



'D2^''S, + ^^S^ 



Ml\a) 

Ml\a) 
Ml\a) 



1 



1 



27ra, 



2( 

m + ra-i m + m2 
2(mi + m2) . , r 



7ni{m + m2) 



3(^ + . ' 



)A2- 






m 


'.(^ 


1 


m 


mi 


_^ 


-B33 


m 



m + rrii m + m2 



mi(m-|-m2) ^/3 mi m 



^D2^^S, + ^So 



Ml\a) = —Ml\a;'D2^'S,+'S,) 



5V3r^ 



2m(3mi-m2) + 2mi(mi+m2) . , /kA . ^ ^ n 
+ + + ^/2(- + — )5i 



m mi 



Mf(a) = Mf (a;i2)2^'5i+35i) 
v2 



Ml\a) 



2'Kas 
Sm^ 



m(3mi-m2) + mi(mi+m2)^ ( + ^ )B 

mi{m + nil) {m + 1712) ^ \/3 ™ 
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Ml\a) 
Ml\a) 

Mf(a) 
Ml\a) 



V2 



Ml\a]^D2^'''S,+''Si) 



1 



2TO(mi-3m2)-2mi(mi+m2) . , 3 1 
■ ■ A2 + Vi( jiill 

mi{m, + nil ){in + 1712) nii m 



mi (mi + ma) - m(mi - Sma) ,^^3/^ L)_B 

mi(TO-|-TOi)(m-|-TO2) ^ mi 3m 



5m2 



70 
27ra 



35M2 



m(mi - Smz) - Wi (wi + ma) -v^fJ^ —)B 

mi{m + nil) {m + 1712) ^ mi 3m 



Mf(a) 



'HMf(a;^7?a^^5i+^5i) 



Mf(a) = -2- 



27rQ!. 



mi(m + ma)^^ + 7lt"m)^^^. 



• ^Ds^'So + 'Si 



M!\a) 

Ml\a) 

A.4 + 5 



^^Ml\a;'D2^'Si + 'Si) 



2'Ka, 



■ ; A2 + V3( ^)-B33 

mi(m-|-mi)(m-|-m2) mi 3m 



M^'{a) = M^^'{a) = 



/-m(mi+m2 + 2m) . , „ ' 

V 6- ; ' -^1 -|-i>01 

(m + mij(m + m2) 



Mf = - 



2nas 
15M2 



2V35mM4 



mi(mi-m)+m2(mi-|-3m) /t^ . , , I x ^ 

;^ ; ; V0^1 + (, 1= J-D()l 

mi(m-|-mi)(m-|-m2) m V2mi 

V2m{mi+m2 + 2m) , , 
A — Jd 



{m + nil) {m + 1712) 



= 



Mt 



3\/35M2 



67r6 



10\/7mM4 



67ra. 



15\/7M2 



\/2[mi (mi - m) + ma (mi + 3m)] _^ _3 ]_ , 

mi(m-|-mi)(m-|-m2) mi m 

2\/7m(mi + ma + 2m) 
(m + mi ) (m + m2 ) 

2\/7[mi(mi-m)-|-m2(mi + 3m.)] 3 1 

-^3 +3( )iJ43 

mi m 



^3 -3a 



mi (m + mi ) (m + m2 ) 
where 

A' =pfPBmy^PAc{lA) + 2pAciO,0)] 

B' =PB(0)[4x/3t;Ac(2, 1, -1) + V6t;^c(l, 1,0) - \/6^;^c(l,0,-l) - 2\/2t;^c(0, 1, 1) + 2\/2t;xc(0,0,0)] 
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A.5 F^S + S 



21M4 



-( 



1 



m + mi 



1 ^ . 1 

— )A3 + — B44 

m + 1712 ^Tn 



3V21M2 



2'm{mi — 3TO2) — 2mi(mi +TO2) 
mi (m + TO 1 ) (to + 771-2 ) 



-3(— --^)i344 

TOi OTO 
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